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Abstract 

Several complete systems of integrability conditions on a spin chain Hamiltonian 
density matrix are presented. The corresponding formulas for i?-matrices are also given. 
The latter is expressed via the local Hamiltonian density in the form similar to spin 
one half XXX and XX Z models. The result is applied to the problem of integrability 
of SU(2) x SU(2)- and SU(2) x [/(l)-invariant spin-orbital chains (the Kugel-Homskii- 
Inagaki model). The eight new integrable cases are found. One of them corresponds to 
the Temper ley-Lieb algebra, the others three to the algebra associated with the XXX, 
XX Z and graded XX Z models. The last two .R-matrices are also presented. 

1 Hamiltonian density and i?-matrix 

It is well known [I],[3], that an i?-matrix, which satisfy the Yang-Baxter equation (in the braid 
group form), 

R 12 (\ - [x)R 23 (\)R 12 (fi) = R 23 (fx)R 12 (\)R 23 (\ - //), (1) 
(usually denoted by R(X)) and the initial condition, 

R(0) oc I, (2) 

where I is an identity matrix, produces according to the formula: 

H=^R(X)\ X=0 , (3) 

the corresponding to an integrable spin chain Hamiltonian density. Let us remain that both 
H, and -R(A) are M 2 x M 2 (M = 2, 3, 4, ...) matrices. For such a matrix X the corresponding 
M 3 x M 3 matrices Xi 2 and X23 are defined as follows: 

X 12 = X <g> I M , X 23 = I M ® X. (4) 



Here Im is the identity M x M matrix. The Hamiltonian may be obtained according to the 
formula: 

H = J2 H n,n+l- (5) 

n 

A matrix R(X) is important for applications in statistical physics. However for the problems 
related to spin chains just the matrix if is a starting point. By this reason it is necessary 
to be able to reverse the formula (3) and get the .R-matrix from a Hamiltonian density. The 
system of necessary conditions for the such operation is an infinite sequence of integrability 
conditions. The first one (the Reshetikhin condition) 

[H l2 + H 23 , [H 12 , H 23 \] = K 23 - K 12 , (6) 

is known long ago [3], the next two are presented in |5j,[6j, while the others are not yet 
obtained. Let us notice that the condition (6) will be satisfied when it will be found an 
appropriate matrix K for the right side. 

In [6] the system of integrability conditions was studied in the context of a power expansion 
of a .R-matrix (l)-(3) with respect to A. A new method for solving the Yang-Baxter equa- 
tion based on an analysis of such expansion was also presented in this paper. This method 
turned out to be rather effective. In particular it allowed to find all .R-matrices for integrable 
spin-ladders [7J. At the same time taking into account only the necessary integrability con- 
ditions this method does not guarantee an existance of a solution. Consiquent solution of 
the integrability conditions together with expanding of quotents for the .R-matrix entries into 
power series gives a possibility at a some stage to suggest their explicit forms. The latter may 
be expressed in the framework of rational, trigonometric (exponential) or elliptic functions. 
Nevertheless the Eq. (1) must be checked for the obtained .R-matrix. 

From the other hand it is better to the reverse the formula (3) having in a stock of knowledge 
a rather big number of sufficient integrability conditions, related to explicit formulas expressing 
an .R-matrix from the corresponding Hamiltonian density. As an example we may suggest the 
sufficient condition virtually obtained in [8] and related to the Temperley-Lieb algebras. [9]. 
Namely if the matrix H satisfy the following system conditions: 

H 2 = —dH, H 12 H 23 Hi 2 = H 12 , H 23 H i2 H 23 = H 23 , (7) 

then the Eq. (6) turns into identity for K = 2H. Moerover the following i?-matrix, 

R TL (X) = Jsh(A + ??) + HshX, (8) 
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where chrj = d/2 according to (3) generates a Hamiltonian density, that differs from H by the 
item Ichrj, which is proportional to the identity matrix and results to an unsufficient shift on 
a constant in the spectrum of the Hamiltonian. 

In the next sections we shall represent the two more sets of sufficient conditions on the 
Hamiltonian density. These conditions were obtained by studying the structures of .R-matrices 
for XXX, XX Z and graded XXZ models of half spin. 

2 The algebras related to Hamiltonian densities 

The following system of relations (satisfied for the XXX model), 

H 2 =oc J, H 12 H 2 3Hi2 = H23H12H23, (9) 

is related to the i?-matrix: 

R XXX (\)= V I + \H. (10) 

The i?-matrices for XX Z and graded XX Z models are expressed from the related Hamil- 
tonian densities according to the formula: 

R XXZ (X) = HshX + ^(H 2 -I) + -^(chV - H 2 ). (11) 
sh?7 shr/ 

In this case the matrix H depends on the real parameter a = chrj (the case a < 1 corresponds 
to the sin — cos R- matrix [3j), and satisfy a system of relations which generalises of (9). First 
of all it includes the Yang-Baxter equation for the Hamiltonian density, 

H 12 H 23 Hi2 = H23H12H23, (12) 

the relation 

(H l2 - H23) 3 = (a 2 + 2)(H 12 - H 23 ), (13) 
and the following set of relations, 

Hi2H 2 3(H 2 2 - I) = (#23 _ -0 H12H23, 

H 12 (H 2 3 - I)H 12 = H 23 (H 2 2 - I)H 23 , 

(#12 — I) H23H12 = H 2 3Hi2(H 2 3 — I) , 

H 12 (Hl 3 -I)(Hl 2 -I) = (Hi 3 -I)(Hf 2 -I)H 23 , 

(H 2 12 -I)H 23 (H 2 12 -I) = (H 2 23 -I)H 12 (H 2 23 -I), 

(H 2 2 -I)(H 2 23 -I)H 12 = H23(H 2 12 -I)(H 2 23 -I), 

(H 2 2 -I)(H 2 23 -I)(H 2 2 -I) = (H 2 23 - 1)(H 2 2 -I)(H 2 23 -I), (14) 
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appearing from a direct substitution of (11) into (1). 

Evidently the system (14) is overfulled and we have written it only to make clear the proof 
of the Eq. (1) for the i?-matrix (11). However it is not so obvious but in this case the Eq. 
(13) is necessary as well as (12) and (14). The following two conditions are also satisfied for 
the standard S —1/2 XX Z i?-matrix: 

H 12 H 2Z = #23#12; ( 15 ) 

H 3 = aH 2 + H-aL, (16) 

Moreover the system (14) may be easily proved from (12), (13), (15), (16). For the i?-matrix 
of the graded XXZ model [3] the relations (12)-(15) are still satisfied, but the Eq. (16) 
substitutes by a more general one: 

H 4 = (a 2 + 1)H 2 - a 2 I. (17) 

In this case an appropriate choice of the minimal set of relations on a Hamiltonian density 
still remains an open question. 

3 Hamiltonian densities and related i?-matrices for spin- 
orbital chains 

The following Hamiltonian density was suggested in |10|-[T2~] in order to describe magnetic 
systems with an additional orbital degeneracy in electron systems of atoms, 

H n , n+1 = (s n s n+1 )( ai + a 2 {r*T* +1 + t*t% +1 ) + a 3 T*r z H+1 + ^a 6 (r* + r* +1 )) + 

+ + t*t% +1 ) + a 5 r*< + i + 2°r(^ + (18) 

Here the spin and pseudospin operators s n and r n are related to the Pauli matrices, 

s k = \ k ®I 2 , T k = ^I 2 ®o- k , k = x,y,z. (19) 

In below the set of coefficients {a 1; ..a 7 } will be denoted by S. The martrix H e End((C 2 ® 
C 2 ) ® (C 2 <8> C 2 )) related to (18) is invariant under the following symmetry, 

H ^Q® QHQ' 1 ® Q' 1 , (20) 

where Q G SU(2) x U(l) and the group U(l) is generated by the matrix r z . 
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The following substitution of coefficients: 

{ai, a 2 , as, a 4 , as, a 6 , aj} — > {ax, a 3 , a^, as, — a 6 , —07}, (2f ) 

relates to the transformation (20) with Q = l2®o x and therefore does not change the spectrum 
of the Hamiltonian as well as does not fail the integrability. 
The similar is true for the following substitution: 

{ai, a.2, 0,3, 04, 05, <2q, aj} — > {ax, —a2, as, —a^, a§, a§, aj}, (22) 

related to the symmetry (20) 

H2n,2n+1 ~~ > ^T2rflln,2n+\ T 2m ^2n-l,2n ~~ > ^ T 2n^'in-\,2n T 2n- (^3) 

for a chain with even number of sites. 

The only one integrable case of the model (18) was known up to now. It relates to Sq = 
{1,4,4,1,1,0,0} p3]-[T5]. In this, S'f/(4)-case the shifted on 1/4/ Hamiltonian density is 
equal to the permutation operator in C 4 ® C 4 , so all the conditions (9) are satisfied, and 
the i?-matrix relates to the forlmula (10). Let us notice that the dual model obtained by 
the transformation (22) is also integrable with the i?-matrix related to the same formula. A 
modern discussion of the S"?7(4)-symmetric model is given in [16], [17]. The SU(2) x SU(2)- 
and SU(2) x [/(l)-symmetric models were also discussed in [18j,[19j. 

After substitution of (18) into the condition (6) we have obtained 8 new solutions, (e, ex, £2 = 

±1), 

Sx = {1,-45,-4,5,1,0,0}, (24) 

5 2 = {0,4,0,l,0,4e,e}, (25) 

5 3 = {l,8ei,4,2ei,3,4e 2 ,-e 2 }, (26) 

5 4 = {l,8exA,2ex,-lAe 2 ,3e 2 }, (27) 

5 5 = {0,4, 0,1, a, 0,0}, (28) 

5 6 = {0,4, 0,1, 0,0, a}, (29) 

5 7 = {l,4e,8,2e,l,0,0}, (30) 

5 8 = {2,45,-8,-5,5,0,0}. (31) 

In the case Sx the matrix H — 1/4/ satisfy the Temperley-Lieb condition (7) (d = 4). and 
has the following form [8]: 

H -i/ = b®b, (32) 
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(or in component notation Hij^i = hjbki)- Bivectors b and b, written in (32) as a column and 
row, is more convenient to represent as matrices, 
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satisfying the following condition, 
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(33) 



-el. 



(34) 



The i?-matrix related to Si has the form (8). For e = 1 this model is SU(2) x SU(2)- 
symmetric. As it was mentioned in [IS] , the corresponding ferromagnetic model (with the 
Hamiltonian density —H) corresponds to the critical point, separating 4 different phases. The 
ground state of this model is highly degenerate. 

In the case S2 the matrix H satisfies the condition (9), and therefore the i?-matrix is given 
by the formula (10). In the case S3 the shifted matrix H = H + |J satisfy (9) and the related 
.R-matrix is also repersented by (10). In the case S4 the shifted matrix H = H + \l satisfies 
(9) and its .R-matrix also relates to (10). 

In the case S5 the shifted matrix H = H + jl satisfies (12)-(16) with a = § and therefore 
the related .R-matrix is given by (11). In the case S§ the Hamiltonian density satisfies (12)-(15), 
(17) and the related .R-matrix also has the form (11). 

Finitely we represent the R- matrices, related to S 7 and S 8 , 
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/i = (l + A)(l + 3A), / 2 = 1 + A, / 3 = 1 + 2A, ff i = 3A(l + A), £ 2 = A(2 + 3A) 
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A dependence of these matrices on H is not yet found. 



4 Conclusions 

The two sets of sufficient conditions of integrability and related explicit formulas experssing 
the i?-matrices from the corresponding Hamiltonian densities are presented. These conditions 
were obtained from an analysis of the i?-matrices and related Hamiltonian densities for the 
XXX and XX Z models of half spin. The suggested approach develops an idea virtually used 
by P. P. Kulish for the models related to Temperley-Lieb algebras. The new method was 
applied for obtaining of the new SU(2) x S'L r (2)-invariant and 7 new SU(2) x £/(l)-invariant 
i?-matrices for spin-orbital chains. 

The author is very grateful to P. P. Kulish for the stimulating discussion and helpful 
comments. 
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